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ABSTRACT 



A new asymptotic method is presented for the analysis of the traveling waves in the one- 
dimensional reaction-diffusion system with the diffusion with a finite velocity and Kolmogorov- 
Petrovskii-Piskunov kinetics. The analysis makes use of the path-integral approach, scaling 
procedure and the singular perturbation techniques involving the large deviations theory 
for the Poisson random walk. The exact formula for the position and speed of reaction 
front is derived. It is found that the reaction front dynamics is formally associated with the 
relativistic Hamiltonian/Lagrangian mechanics. 

PACS numbers: 05.70.Ln, 82.20.Db 
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Phenomena of the wave propagation in the non-equilibrium media described by reaction- 
diffusion equations have attracted considerable interest in a wide variety of scientific fields 
including physics, chemistry, biology, etc. The excellent reviews of the work in this area can 
be found in the books [1-4]. Of fundamental interest is the rate at which the wave propagates 
through the non-linear dissipative system. The basic common feature in many examples is 
that the transport process determining the propagation rate is described by the conventional 
diffusion (Fick's law). In this case the propagation velocity u can be found from a simple 



dimensional analysis, that is u ~ vDU, where D is the diffusion coefficient and U is the 
characteristic reaction rate constant. 

It is well known that the diffusion approximation gives rise to the infinite speed of 
heat /mass propagation that is if a sudden change of temperature/concentration takes place 
somewhere in the space, it will be felt immediately everywhere with the exponentially small 
amplitude. It is therefore desirable to have a theory for the nonlinear wave propagation in 
which the boundness of the transport process would be taken into account. The purpose 
of this Letter is to present such a theory giving a new asymptotic method for calculating 
the propagation speed for the traveling wave in the reaction-diffusion system involving the 
diffusion with a finite velocity [5-10] and the chemical kinetics of the Kolmogorov-Petrovskii- 
Piskunov (KPP) type [1-4, 11-17]. 

Our starting point is a phenomenological system of the one-dimensional equations for 
the time evolution of the scalar field p and its flux J 

where U is the reaction rate constant corresponding to the KPP-kinetics, D is the diffusion 
coefficient corresponding to the Fick's law and r is the relaxation time. When U — 0, the 
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system (1),(2) reduces to the telegraph equation [5-10] 



d 2 p dp „<9 2 p . . 

when r = 0, we have the classical KPP-equation [1-4] 

f^g + £M1 _ p , (4) 

If we now solve the equation (2) with the initial condition J (0, x) = we may eliminate 
J from (1) to obtain the single equation for p 

dp D f l ( t - s\ d 2 p(s,x) TT 

This equation may be considered as a generalization of the KPP-equation (4) to the case in 
which the finite speed of the transport process is taken into account (r ^ 0). 
We specify the following initial condition 

p(0,x) = 6(x), (6) 

where 9 (x) is a Heaviside function 9 (x) = 1 for x < and 9 (x) = for x > 0. 

The basic problem is to find the traveling wave solution tp (x — ut) to the problem (5), (6), 
where ip (z) is a monotonically decreasing function such that ip ( — °°) — 1 an d "0 (oo) = 0, 
and u is the speed at which the wave profile ip moves in the positive x-direction. For the 
KPP-equation (4) (r = 0) with the initial condition (6) the traveling wave moves with the 
velocity u = \^ADU [1-4]. We expect that for (5), (6) the speed u = y/lDU f (tU) , where 
/ (z) is the dimensionless function such that / (0) = 1. It should be noted that our method 
of calculation will be nonperturbative in a sense that we do not treat the relaxation time r 
as a small parameter. 

We are interested in the long-time large-distance behavior of the traveling wave solution 
of (5), (6) as t — > oo and x — > oo . It is convenient therefore to make the scaling [12-17] 



and rewrite the Cauchy problem (5), (6) for p £ (t, x) — p (|, | ) in the following form 
5p £ D ft ( t-s\ d 2 p £ (s,x) , [/ F/ FX . „, . 



It is clear from this equation that the scaling (7) describing a simultaneous contraction of 
time and space alternatively corresponds to the rapid chemical reaction and slow transport 
process. We expect that after rescaling the wave profile develops into the reaction front: 
p £ (i, x) — ip {^-f^j tends to a unit step function 9 (x — ut) as e — > 0. 

Our goal is now to find a function G (t, x) determining the position of the reaction front 
that is 

hm p £ (tx) = l° ifG (^) <0 ( 9) 
e^o ' I 1 otherwise. 

In this Letter we restrict ourselves in finding the upper bound for p £ (t, x) in the form 

p £ (t,x) <exp| G ^^ | as e^O. (10) 

It is clear from (10) that p £ (t,x) — > if G (t, x) < and e — > 0. One can prove that 
p £ (t, x) -> 1 if G (t, x) > and e -> 0. 

It follows from the property of the KPP-kinetics in (8) that 

(f(t,x)<<ff{t,x)exv^j), (11) 

where </? e (t, x) is a solution of the linear initial problem 

dif £ D ft _ / (U + r- 1 )(t- s)\ d 2 cp £ (s,x) 



dt 



t/o^H r^H^ s ' = «(*)■ (12) 



Our strategy to find the function G (t, x) is to analyze the above Cauchy problem in 
terms of the probability theory and thereby to obtain an estimate of ip £ (t, x) in the limit 
e — > 0. The basic idea is that we can deal with the problem (12) in terms of the random 
walks of Poisson type [7-9] . If we introduce the notations c (velocity) and v (frequency) such 
that 

9 D 1 

c 2 = — , 2u = U + -, (13) 

T T 



than the solution of linear initial value problem (12) can be written as an expectation value 
of the initial distribution 9 [7,8] 

<f(t,x)=E9(x(t)), (14) 

where E is the expectation operator and x (t) is a random Poisson walk, i.e. a solution of 
the stochastic differential equation 

v(-), x(0)=x, 0<s<t, (15) 



ds \e 

where v (s) is the Markovian dichotomous velocity taking only two values {c, — c} with the 
frequency v [18]. From the probabilistic point of view the key factor underlying the nonlocal 
character of (12) is that the dynamics of x (t) is non-Markovian [18]. To obtain an estimate 
of ip e (t, x) as e — > we need to know an explicit expression for tp e (t, x) as a path-integral 
[19] 

<ff(t,x) = Je(x(t))P[x(.)]Vx, (16) 
where P [x (•)] is a conditional probability density functional for the random process x (s) 



P [*(•)] 



dx / s 
ds \e 



P[v(-)]Vv, 



where 5 [•] is the 5— functional that is the extension of the ordinary 5— function to the func- 
tional integration [19]. By using (16) and the formalism based on the auxiliary function u 
[20,21] we can write down the following expression for ip £ (t, x) 

<p e (t,x) = J J Je (x (t)) exp J* u (s) (^~v j ds^j P [v (•)] VvVuVx. (17) 

In the "weak noise limit" e — > one can get the following estimate for ip £ (t,x) (the 
details of calculation involving the large deviations theory for the Poisson random walk will 
appear elsewhere) 

(p £ (t, x) ~ exp I — - min / f p^- — H (p) ) ds 1 , (18) 
v ' { e x(o)=x-,x(t)=o Jo \ ds v 7 J v ' 
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where the function H (p) has the form of the relativistic Hamiltonian 



H (p) = c\J m?c 2 + p 2 — v (19) 

with the "effective mass" m = vc~ 2 . One can rewrite (18) in terms of the relativistic La- 
grangian [22] 

ip £ (t, x) ~ exp | — - min f Lds 1 , (20) 



E x(0)=x;x(t)=0 JO 



where 



L = —mc 2 



1 ( dx x 2 



f"-\Ts) + "' < 21) 

We are now in a position to complete the derivation of the function G (t, x) determining 
the reaction front position and its speed. One finds after straightforward calculation that 
the optimal trajectory giving the minimum in (20) is x (s) = —j.s + x and the corresponding 



minimal action is mc 2 t\j 1 — ^ (j) . By using (10), (11) and (20) and the relation m = vc 
we obtain 



G (t, x) = Ut - vt + 1 - — (^j . (22) 
Equating G (t, x) to we obtain the position of reaction front x (t) 



v 



x(t)=ut, u = c\l-[ , U<v. (23) 



Taking into account (13) the speed of the reaction front u can be rewritten in terms of the 
phenomenological parameters D and r 



u 



rU<l. (24) 



\(1 + tU) 



It follows from (23) and (24) that the speed u takes the maximum value c = y^, when 
tU = 1 or U = v. If the velocity of propagation c were infinitely great and the time r were 
infinitely small such that D = c 2 t = const (the diffusion approximation for the random walk 
of Poisson type), equations (23) and (24) would merely give u = y/ADU - the classical result 
of the KPP-theory [1-4]. 
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In summury, we have derived the exact formula for the position and speed of the reaction 
front in the one-dimensional dissipative system involving the diffusion with a finite velocity 
and the KPP-kinetics. It has been found that the reaction front dynamics for such a system 
can be formally associated with the relativistic Hamiltonian/Lagrangian mechanics. There 
are several possible directions to explore by the method developed here. First one may 
study the influence of non-uniform distribution of the reaction rate constant U which might 
induce the jumps of reaction fronts [11]. One can also extend the analysis to describe the 
interaction between the turbulent diffusion with a finite velocity [10] and the KPP-kinetics 
in the three-dimensional space [17]. 

The author would like to thank Rupert Klein and Alexander S. Mikhailov for interesting 
and helpful discussions. The research was supported in part by DFG Project KL 611/5-1. 
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